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We show that even a rather minimal extension of the Einstein-Hilbert action by a nonminimal 
coupling of the scalar field to the Ricci curvature scalar results in configurations that resemble more 
the dark energy stars then the ordinary boson stars. Even though many of those configurations 
are endowed by negative principal pressures, the strong energy condition, as a signal of repulsive 
gravity, is not significantly violated in these configurations. When imposing restrictions on matter 
1 from energy conditions we find that the maximally allowed masses are shifted to the lower values 

y—\ ' due to the violation of the weak and dominant energy conditions. We also calculate the effective 

compactness and show that its maximum value is attained in the region of negative pressures, and 
is greater then that in ordinary boson stars. Moreover, we develop a universality technique which 
allows to efficiently map small configurations, that are easily solved by numerical methods, to large 
astrophysical objects. 

Q: 

\Q '. I. INTRODUCTION 



i— i ■ One of the most peculiar predictions of classical general relativity, al least if matter obeys the strong energy 
condition (SEC) are black holes. (According to the SEC the sum of the energy density and pressures, p+J2 Pi — 0' 
cannot be negative.) When quantum effects are included, black holes lead to a number of thermodynamic 
paradoxes associated with Hawking radiation and the implied information loss in black hole spacetimes, thus 

, questioning whether the final stage of a massive star collapse is a black hole, or perhaps some other as-yet- 

not- understood dense object, that stops a further collapse. Sakharov was the first that introduced the concept 

t-H | of nonsingular collapse through the equation of state for the cosmological dark energy (for which p ~ —p) as a 
super-dense fluid [l| and then Gliner assumed that such a fluid could be the final state of gravitational collapse Q • 
Inspired by these ideas Mazur and Mottola investigated alternative configurations which led to a solution dubbed 
gravastar (gravitational vacuum star) 3]. This anisotropic, highly compact astrophysical object consists of a de 
' Sitter core and through vacuum phase transition layer matches an exterior Schwarzschild spacetime by avoiding 
an event horizon formation. Although gravastar configurations rest upon a very attractive idea, all these models 
are macroscopic in the sense that their foundation rest on studying Einstein's theory in presence of a matter 
fluid that obeys some phenomenological equation of state, and as such do not have a proper field theoretic 
foundation. Both cases - when the energy density is distributed on thin-shells Q and when it continuously varies 
• • ' throughout the star 

043 - rely on the so-called ^4nsd'tee-approach. In this approach Einstein's equations are 
solved in presence of a radially distributed matter fluid, for which an equation of state or some other relation 
among the thermodynamic functions (the energy density, the radial and tangential pressures) is provided. All 
these models are essentially toy models, and they are important in the sense that they can be used to provide a 
better understanding of the main characteristics of black hole mimickers. But, a complete understanding of these 
objects will be attained only if we can provide faithful microscopic (field-theoretic) models for these objects. 

Apart from a better understanding at the fundamental level, field theoretic models of highly compact non- 
singular objects, obtained from a suitable lagrangian of interacting matter fields, can provide a fundamental 
explanation for the anisotropy in the principal pressures, which occurs naturally in the stars made of scalar fields, 
the so-called boson stars. Boson stars are nonsingular asymptotically flat solutions of the Einstein-Klein-Gordon 
field equations which govern massive complex scalar fields coupled to gravity. The extensive research started by 
Kaup [8J, who has introduced the notion of the gravitationally bound state of scalar particles. Soon many papers 
considering various versions of scalar field configurations appeared |9r|l4|. The growing importance of boson 
stars resulted in extensive research which has been reviewed in |15l4l8l |. When the boson star configurations 
are considered, one immediately recognizes that a massive scalar field, even if self-interacting, cannot produce 
anisotropy which could support an object with (asymptotically) de Sitter interior. Albeit boson stars belong 
to the realm of very compact objects, it turns out that getting closer to the main black hole features requires 
modification of general relativity. Even though Einstein's theory has passed all observational tests in the weak 
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field limit, the true theory of gravity may differ significantly in the regime of strong gravitational fields. Moreover, 
large scale cosmological observations and conceptual difficulties in quantizing general relativity call as well for its 
modifications. 

In this paper we show that even a rather minimal extension of the Einstein-Hilbert action by a nonminimal 
coupling of the scalar field to the Ricci curvature scalar results in configurations that resemble more the dark 
energy stars then the ordinary boson stars. Even though many of those configurations are endowed by negative 
principal pressures, the strong energy condition, as a signal of repulsive gravity, is not significantly violated in 
these configurations. Yet, the maximum effective compactness is attained in the region of negative pressures, and 
is greater then that in ordinary boson stars. This fact supports an idea that the dark energy stars might present 
a promising black hole mimicker. While some attempts have been made to study dark energy stars (see for 
example Refs. [Hrllij]), which are loosely speaking in literature taken as objects that contain a negative pressure 
somewhere in the interior, no systematic study has been so far performed of whether and when boson stars in 
nonminimal setting violate energy conditions. In this work we fill that gap. 

The paper is organized as follows. In Sec. II we present the basic Einstein equations for spherically symmetric 
configurations of a nonminimally coupled complex scalar field. In Sec. Ill a brief description of ordinary boson 
star solutions is presented which leads to a situation in which a necessity of additions mechanism is needed. The 
nonminimal coupling is introduced in Sec. IV, where we present solutions with required anisotropic behaviour of 
pressures, i.e. dark energy- like stars comprising negative pressures. In that section we also perform analysis of 
the parameter space for which the weak and dominant energy conditions are violated. Moreover, we investigate 
the effective compactness. Finally, in Sec. V we discuss our results and give an outlook for future work. 



II. THEORY BEHIND - EQUILIBRIUM CONFIGURATIONS 

For gravity we take the standard Einstein-Hilbert action: 



Seh = d xy^g 



R 



16itGn 



(1) 



where Gjq is the Newton constant, R is the Ricci scalar, and g is the determinant of the metric tensor g^ v , which 
is given by 



g^=dmg{-e^ r \e x{r \r\r 2 



sin 



(2) 



The space-time metric is static and spherically symmetric as we are interested only in spherically symmetric 
equilibrium configurations. For matter we take an action of a complex scalar field with a mass and a quartic 
self-interaction coupled nonminimally to gravity: 



<\> - m l<l>* "<t> - 



+ ^R<f>* 



(3) 



where £ measures the strength of the coupling between scalar field 4> and gravity via the Ricci scalar R and cf>* 
is the complex conjugate of <j>. It is worth noting here that in order to produce stable configurations, the scalar 
field must be complex. According to the Derrick's theorem (25j regular, static, nontopological, localized scalar 
field solutions cannot be created by real scalar fields (see [26| and e.g. [HI). 

The energy-momentum tensor of a complex scalar field is obtained by varying its action with respect to the 
metric tensor g^ v : 



T% = IS^Sfard^-g^u 

- 2^*4>G fil/ + 2£V Al V t/ (</>» - 2^D(0V 
By varying now the full action 

S = Seh + S$. 

with respect to the metric tensor g^ v we obtain the Einstein equations: 



(4) 
(5) 



(6) 
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The Klein-Gordon equation, the equation of motion for the scalar field 4> ( or 
identities or by varying ([5]) with respect to <jf (or <p), resulting in: 

p-ml-\4<l>*<t> + ZR]<t> = 0. 

In order to proceed we choose a harmonic time-dependence for the scalar field 

<j>{r,t) = <t> {r)e- wt , <Po(r) G K. 



is obtained from Bianchi 
(7) 

(8) 



Even though the scalar field that induces the metric is time-dependent, the energy momentum tensor created by 
this field is time-independent and thus leads to time-independent metric functions. Hence, condition (jHJ does not 
contradict the Birkhoff theorem. Furthermore, the same Ansatz for the classical field was also used in Ref. [27l | 
bearing the name coherent state, presumably alluding to their resemblance to quantum coherent states. 1 Highly 
excited field configurations were used to explain flat rotation curves inside galactic halos in Ref. [281 ] . Furthermore, 
in Ref. [2^ | it was shown that it is possible to construct a stable multistate boson star, with coexisting ground 
and first excited states. 

Upon inserting the Ansatze (JSJ and §2§ into ((4]) one gets for non- vanishing components of the stress energy tensor: 
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Similarly, the scalar field equation of motion (J7J becomes: 

v' — \ r 



00 + (; + V -^~) 0o " e A (ml + X^l - u 2 e~ v - o 



0. 



(9) 
(10) 



(11) 
(12) 

(13) 



By virtue of (IT31 it is possible to eliminate the second derivative of the scalar field O ' in the components of 
the energy-momentum tensor (|9l fT2)) . leading to the following form for the first two Einstein equations [G" — 
8irG N T»): 



[1 + 2£(8 7 rG JV )0o] G\ 



8nG N {( -ml- u 2 e-» - -±<fi} <g - e - A (l + 4£)0 O 

-4£ [ml + X^l - u?er v - £R] <g + 2^ x i/ '0o0 o }, 
[l + 2^G N )ct>l]G r r = 8nG N {(-ml + w 2 e~ v -^fjeft + e^fi 2 - 2&~ x / / ^ 



where 



Gt = 



-A 



Gl = e- x 
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(14) 
(15) 

(16) 
(17) 



One should keep in mind however that, a scalar field written as in {5), apart from the ground state, can also represent excited 
states with higher energy, and a particular combination of these states can indeed form coherent states. In general, these states 
contain 'coherent' radial oscillations, but do not in the usual sense constitute quantum coherent states. 
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There is one more independent equation. Instead of using the (99) Einstein equation (or the equivalent (ipip) 
equation), it is in fact more convenient to use the trace equation, = —R = 8t:GnT^, leading to: 



R = 



8nG N {2mU 2 + 2(1 + 60 le^^ 2 + (m 2 - to 2 e- + A^ 2 )l 2 } 

1 } - (18) 



1 + 2£(1 + 608ttG 



N<P 



It is instructive to add a couple of remarks on this equation. For the case of conformal coupling, £ = —1/6, the 
only non- vanishing term in the Ricci curvature scalar is the scalar field mass. Hence in limit of a vanishing scalar 
field mass, for which the Ricci scalar is zero, one obtains a conformal gravity limit, as expected. Nevertheless, 
for £ =^ —1/6, as we shall see in the subsequent sections, a variety of configurations is possible. 
Equations (fT4HT5| and (fT8|) constitute the central equations in this work. 

A. Dimensionless variables 

Before we proceed to solving Eqs. (114lTl5l) and (TT51) . for the purpose of numerical studies, it is convenient to 
work with dimensionless variables/functions. To this purpose we perform the following rescalings: 

► x, 8irG N 4>o{ r ) 2 -> c(r) 2 , 



N 



8ttG n R -> R, 8nG N m^m^ 8nG N w z -> Q 2 . (19) 
Upon these transformations all variables/functions get expressed in terms of reduced Planck units: 



he 

fhp = top/V8tt, with the Planck mass m P = \j — = 1.2209 x 10 19 GeV/c 2 = 2.17651(13) x 1CT 8 kg, 

Gn 



/ hG 

Ip = VOwlp, with the Planck length l P = J — ^ = 1.616199(97) x 10" 35 m. (20) 
The rescaled (dimensionless) differential equations to be solved are then: 

A' 



1 - e x e x (rhl + tfe~ u + ^a 2 )a 2 + (1 + 40a' 2 - 2&'aa l 



x 1 + 2£ct 2 

4x£e x (rh 2 - Q 2 e~ v + X^cr 2 - £R)a 2 
+ 1 + 2^ ' 

(e A - 1)(1 + 2£<t 2 )/.t + xe x {-m 2 + Q 2 e~ u - \a 2 )u 2 + xa' 2 ~ 8^aa' 

1 + 2£cr 2 + 2£xaa' ; 

2 V — A \ j s 2 2 2 —v 



(21) 
(22) 



^x 2 
with the dimensionless Ricci scalar 



d' + e A (m 2 + A <r 2 - t^e^ - £i?)cr, (23) 



R = 



2m/a 2 + 2(1 + 60 (m 2 - £ 2 e^ + \^a 2 )a 2 + er x o 12 



l + 201 + 60cr 2 



(24) 



where now the primes denote derivatives with respect to x. 

Equations (I21H23[) yield a unique solution (that depends of course on Co) when subject to the boundary conditions: 

(1) A(0) = 0, (2) v(oo) = 0, (3) cr(0) = cr , (4) cr(oo) = 0. (25) 

The first boundary condition ensures that the mass function m(r) defined in terms of the metric function as 

i 2m(r) 2rh(x) 
e~ x = 1 ^ = 1 K —t (26) 
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is zero at r = (or equivalently at x = 0). The second boundary condition in (|25[) ensures asymptotic flatness at 
large distances, 

The third and fourth boundary conditions in (|2~51) are typical for boson stars with a positive scalar mass term 
(m 2 , > 0). 2 Equations (|2T1[23| together with the boundary conditions (f25j) constitute an eigenvalue problem 
for uj - that is, for each central field value Co there is an unique u) that satisfies the given boundary conditions. 
The ground state is characterized by zero nodes in the field a(x) (defined as the points x where <j(x) — 0), while 
the n-th excited state has n-nodes in o~(x). In this paper, if not explicitly stated otherwise, boson stars in their 
ground state will be studied. 

We solve these nonlinear, mutually coupled, differential equations numerically by using the software package 
Colsys [H. 



1. 



(27) 



B. Universality 

In order to solve the problem numerically we need to specify the set of parameters {A^, m 2 , £> 2 , £}. However, 
for a successful numerical integration these parameters cannot be very different from unity. On the other hand, 
in physically interesting situations these parameters may wildly differ from unity. For example, compact stars 
have radial size that is measured in kilometers, while numerical solutions give objects whose size is of the order of 
the Planck length, lp ~ 10~ 38 km, obviously not very useful. In order to overcome this impasse, we observe that 
the dimensionless equations (pHHlMf possess a 'conformal' symmetry. Indeed, Eqs. (f2"Trf2~3| are invariant under 
the following conformal transformations 



A - R 



;,2 



x^(3x, A -> R ^>jp, ™ 2 ^-^-> ^^^2' ( 28 ) 

How the mass of the whole boson star changes due to these rescalings can be estimated from the identity 

M~pR 3 , (29) 

where p is the density which can be approximated by the value of the potential at a scalar field maximum 

p ~ V(a ) ~ m 2 a 2 + A^. (30) 

On the other hand, from the virial theorem, according to which star's gradient energy ~ potential energy, the 
radius of the star, i.e. its core in which most of its energy is contained, can be estimated from 



(v^) 2 ~ v (ft, 

m!a 2 + A a 4 . (31) 



a 2 



R 2 

This then implies that the mass of a boson star scales as the radius, M oc R, leading to 

M -4 (3M. (32) 

For example, for a compact object whose radial size is R ~ 10 km = /3lp, we obtain that (3 is of the order of 
(3 ~ 10 38 . It then follows that the mass of the scalar field changes from ~ mp to ~ 10 _38 mp and the 
coupling constant from ~ 1 to \,p ~ 10 -76 . In light of Eq. ([20|). the total mass from M ~ fhp changes to 
M ~ 0.2 M Q , where M Q = 2 x 10 30 kg is the solar mass. 

On the other hand, one can start by setting the scalar field mass and estimate the resulting star radius and 
its total mass. This allows one to build models that can account for astrophysical obj ects of vastly different sizes, 
namely from dark compact objects (3^, [UJ to galactic dark matter halos [13, [H, l33j . 



2 On the other hand, the appropriate boundary conditions in the case of a negative mass term (m^ < 0) are cr(0) = 0, cr(oo) = co, 
resulting topologically stable configurations known as global strings. More generally, a multicomponent scalar field with appropriate 
boundary conditions can generate global topological defects which have been extensively studied in cosmology. 
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FIG. 1: The star mass as a function of the central field value ao and X^, for £ = on the left panel and the maximal 
effective compactness as a function of on the right panel. Also mi = fh%. 

III. ORDINARY BOSON STARS: CASE OF MINIMAL COUPLING 

Since the properties of the boson stars with quartic self- interaction are quite extensively studied in Ref. [10( , here 
we shall only briefly discuss their main characteristics. Perhaps the most peculiar feature of these configurations 
is the anisotropy in their principal pressures. Whereby in the (usual) fluid approach to compact stars, anisotropy 
is treated as a rather dubious and speculative concept (e.g. in the physics of neutron stars), it appears as a fairly 
natural property of boson stars. One can verify this by inspecting Eqs. (| 1 Oti 1 1 p . which for £ = yield: 

H = Pt - Pr = -2e- X (/) / 2 . (33) 

Here we have identified the components of the energy-momentum tensor as 

T l = diag(- p,p r ,Pt,Pt), (34) 

where p is the energy density, p r the radial pressure and pt is the tangential pressure (pt = pe =p<j>)- 
From Eq. (|33j) we see that anisotropy is strictly a negative function of the radial coordinate. This fact entails 
that, regardless of the coupling strength, for minimal coupling, one can only create configurations with p r > p t . 
In the next section we elaborate more on the consequences of this fact. 

In order to build a viable astrophysical object, its stability is clearly a basic requirement. Stability of boson 
stars has been extensively studied in the literature both analytically [341436T ] and numerically [37l - l39l | 3 . Numerical 
methods include dynamical evolution of the system at hand, whilst the analytical one rely on the standard 
Chandrasekhar's methods, i.e. studying the response to a linear perturbation of static equilibrium configurations, 
whereby the total particle number is conserved 4 . Both avenues, however, lead to the same conclusion that can be 
summarized as follows: there exists a critical value of the central field a c for which the ground state of boson star 
(nodeless in cr(x)) will be marginally stable upon small radial perturbations. For this critical field value, the total 
mass of the star exhibits turnaround in M(<ro)-curve (see e.g. [41]). Then the configurations left from the peak 
are stable and those right from the peak are unstable leading to the collapse to a black hole or a dispersion at 
infinity. An interested reader may find a discussion in Ref. [42j on what is the likely fate of the boson star for the 
right-from-the-peak configurations in the M(tro)-curve. In the absence of the self-interaction it was found that 
the maximally allowed mass is M max — 0.633 m| lanck /m^,, and by switching on the self-interaction it increases 
as M max = O.22v / Ampianck/'Ti0, where A = X^A-kGntti 2 ^). In the left panel of Fig. Q] we show the star mass as 
a function of the central field value exo and X<j>- As the amplitude Co increases its mass also increases (while its 
radius decreases). For increasing A^ the maximum mass also increases while the critical central field value <j c 
decreases. 



3 The catastrophe theory is another interesting method that can be found in Ref. |4C|| . 

4 Since the action is invariant under the U(l) symmetry, according to the Noether theorem, there is a conserved (scalar) charge 
density. 
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FIG. 2: Left panel: the scalar field as a function of the radial coordinate and in the inset the metric functions g t t and 
g rr - Right panel: the energy densities and the principal pressures in the inset. The solid curves are plotted for A^, = 
({cr c , M max } = {0.27, 0.633 m P }) and the dashed curves are plotted for A = 100 ({a c , A/ max } = {0.095, 2.257 rh P }). Also 
m\ = mp and £ = 0. 



Albeit self-interacting boson stars have been extensively studied in the literature, the explicit behaviour of 
their thermodynamic functions, i.e. of the energy density and the principal pressures, is lacking. Thereby, here 
we plot these (and other relevant) functions for two distinctive regimes: a vanishing self- interaction and a large 
self-interaction coupling constant. For any given value of the coupling constant there is only one configuration 
that meets those of the maximally allowed mass. In Fig. [2] we show two such configurations for = and 
= 100. On the left panel, the profiles of the scalar field and the metric functions (inset) are shown, while 
on the right panel we show the behaviour of the energy densities and the corresponding pressures (inset). Two 
main criteria can be read off from these graphs. The first is the interplay among the central field value and the 
radius: while one is increasing, the other one is decreasing and vice versa. An important consequence of this 
trend is equivalence between the M(<jq) and the M(i?)-curve. That is, both curves exhibit turnaround behaviour 
for equal maximally allowed masses, and hence either can be used for stability analysis. Second, the anisotropy 
in the principal pressures becomes less prominent due to the inclusion of self-interaction. This behaviour implies 
that boson stars built from strongly self-interacting fields tend to be more isotropic. Indeed, in this regime, boson 
stars behave like a polytrope with the equation of state pt ~ p r oc p 1+1 / n , where n is the polytropic index. 
It is also worth noting here that the negative anisotropy (1331) cannot rise to more exotic structures with negative 
principal pressures, found in dark energy stars (e.g. gravastars). For the latter, one needs anisotropy to be a 
positive function of the radial coordinate (see e.g. Refs. [H, HH). This is the main reason why we extend this 
analysis to include nonminimal coupling. 



IV. DARK ENERGY-LIKE STARS: EFFECTS OF NONMINIMAL COUPLING 

Spherically symmetric static configurations of a nonminimally coupled scalar field modeled by the action ([3]) 
in the absence of the quartic self-interaction were studied by Bij and Gleiser in Ref. [4^ . Adopting the M(oq) 
stability criterion, the authors calculated the critical (maximally allowed) mass and the critical particle number for 
a variety of values of the coupling constant £. The analysis is performed for boson stars both in the ground state 
(no nodes in the scalar field) as well as in excited states (higher nodes in the scalar field). However, the authors 
did not analyze the behaviour of the thermodynamic functions, namely of the energy density and pressures. 
In the case of nonminimal coupling, the anisotropy becomes rather convoluted function of matter and geometry 

n = -2e- A 0o 2 - 2£(G« - G r r )$ + 2^ x (V + ^ 

- 4£ (ml + X^l - wV - £R) 4>l- (35) 

As mentioned in the previous section, it is likely that Eq. (|35|) may become positive for some radii, which is an 
important ingredient of building microscopic configurations with negative pressures. 

However, when dealing with spherically symmetric, localized, configurations of matter it seems reasonable to 
invoke the energy conditions as important criteria for physically acceptable matter. 
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FIG. 3: The energy density as a function of the radial coordinate. The solid curve stands for £ = = —1, the dashed 

curves are for £ = —2, —4, —6 from left to right, respectively. Also = and = m 2 P . 

A. Constraints from the energy conditions 

Various energy conditions have been proposed as reasonable physical restrictions on matter fields. They origi- 
nate from the Raychaudhury equation together with the requirement that gravity should be attractive (see e.g. 
Refs. [23|, HH, H|| ) . When translated to the energy momentum tensor for an anisotropic matter they read 

The Weak Energy Condition (WEC) p > 0, p + p r > 0, p + p t > 0, 
The Dominant Energy Condition (DEC): p — p r > 0, p — p t > 0, 

The Strong Energy Condition (SEC): p +p r + 2p t > 0. (36) 

The weak energy condition imposes the requirement of a positive energy density measured by any observer. 
Also the energy density plus pressures in any direction needs to be positive. The dominant energy condition 
requires that the pressures of the fluid do not exceed the energy density, so that the local sound speed in any 
observable fluid is always less then the speed of light in vacuum. The strong energy condition has very interesting 
implications. Its violation lead to regions of repulsive gravity such as in cosmological inflation and gravastars. 
Hence it is reasonable to require that the WEC and DEC are satisfied by a fluid, but that the SEC may be 
violated. 

With or without self-interaction it turns out that the WEC is obeyed for nonminimal couplings only if greater 
than a certain (negative) critical value £ > C^f , whereby £^f c decreases very slowly as increases. As an 
example of the indicated transition, we plot the energy density on Fig |3] where it is shown that violation of the 
WEC is more prominent as the value of the nonminimal coupling decreases. 

One example of a spacetime that violates the WEC is that of a wormhole (see e.g. Ref. |47|). Some other examples 
would include a more exotic matter. Although this energy condition is also violated by certain quantum fields, a 
positive energy density is an essential feature of the classical forms of matter. A consequence of the requirement 
that the WEC is satisfied is a shift in the "maximally" allowed masses to lower values as depicted on Fig. |4]for 

f s tWEC 
S ^ Scrit • 

In addition, we also require that the energy is not transported faster than light, and hence the dominant energy 
condition should be satisfied. Another constraint on parameter space emerges from the requirement £ > £®f t c 
as shown in Fig. 2J The dashed curves represent the maximally allowed masses for the by-the-WEC-and-DEC 
modified configurations, while the solid curves correspond to the old (non-modified) configurations. As such Fig. 2] 
represents an important result of this paper due to the fact that it establishes new configurations for stars that 
satisfy the WEC and DEC. 

The strong energy condition is also violated for certain nonminimal couplings, £ > £f r f t C - As opposed to the 
problem of violating the WEC and DEC, a violation of the SEC is actually favorable in building highly compact 
objects. Namely, a region of a compact object that violates the SEC exhibits repulsive gravity, which is desirable. 
Violation of the SEC plays an important role in the early universe cosmology, where it is used to explain the origin 
of Universe's large scale structure generated through matter and gravitational perturbations amplified during a 
hypothetical inflationary epoch in which the SEC is violated. It is also an essential component of gravastars, 
which in their interior, where p r (0) = pt(0) — —p(0), strongly violate the SEC. In the case of gravastars, violation 
of the SEC is crucial for large values of compactness. Unfortunately, here the SEC is significantly violated only 
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FIG. 4: The maximum mass as a function of the coupling £ for A^ = {0, 20, 50, 100} from bottom to top. Also m^, = nip. 
For £ < the dashed curves describe configurations that obey the weak energy condition and for £ > configurations 
that obey the dominant energy condition. The solid curves describe configurations that are not constrained by energy 
conditions. 
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FIG. 5: The energy density and the principal pressures (insets) for m| = mp, £ = —4 and for a) {A^, a c } = {0, 0.050} on 
the left panel and b) {\^,,a c } — {100,0.034} on the right panel. 

if the DEC is violated. Nevertheless we shall explore some effects of violating the SEC in the next subsection. 

B. Energy density and pressures profiles 

It is now of interest to explore thermodynamic functions, namely the energy density and the principal pressures. 
Depending on the strength of the self-interaction, configurations with negative principal pressures emerge, that 
can be described by the equation of state p r oc — p@ . This particular equation of state (EoS) is used to describe 
dark energy stars. Even though these configurations exhibit negative principal pressures, the strong energy 
condition is not violated thus excluding regions with repulsive gravity. 

One such configuration is shown on the left panel of Fig. O As a matter of fact, in the absence of self-interaction, 
for £ < £5Jft C au configurations lying on the M max (£)-curve can be described by the EoS of a dark energy star, 
p r oc — pr . When the self-interaction increases, the pressures increase as well, as can be seen on the right panel 
of Fig. \5\ Nevertheless, no matter how large the self-interaction is, the dark energy star-like configurations are 
obtained by choosing an appropriate (i.e. negative enough) nonminimal coupling. This effect is clearly shown by 
comparing Fig. [B] with Fig. [5j It is also of interest to observe that the transversal pressures of these configurations, 
as positive near surface, are like those of gravastars. This fact brings us to the idea that the gravastars, as not 
yet formulated within the field theories and as such still of interest to explore, might be produced in modified 
gravity that includes higher order terms in the Ricci scalar, Ricci tensor and/or Riemann curvature tensor. 

However, proper dark energy stars, i.e. with negative pressures and violating SEC, can be obtained for A^ < 0. 
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FIG. 6: The energy density and the principal pressures (inset) for £ = —8, = 100, a c = 0.02. Also mi = m 2 P . 
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FIG. 7: The energy density (solid curve), radial (short-dashed curve) and transversal pressure (long-dashed curve) for 
£ = -0.9, A = -30, cto = 0.26. Also m| = trip. 

We present one such solution in Fig. [7] It is interesting that this solution violates only the strong energy condition 
while the weak and dominant energy conditions are obeyed. Nevertheless, theories with negative potentials yield 
Hamiltonians that are unbounded from below, and are at best quasi-stable, i.e. field configurations will eventually 
'decay' into large fields and roll down to infinity, where energy is minus infinity (see e.g. |28|). 

When excited states of these configurations are considered, the energy density and pressures oscillate in space. 
Both pressures are now positive functions of coordinate in the region near the surface thus resembling gravastars 
solutions. We show one such configuration in Fig. El However, even though stability might not be questionable in 
this setting, the weak and dominant energy conditions are violated. Yet, it was argued in Ref. [28| that a galactic 
halo consisting of highly excited states of ordinary boson stars could explain the rotation of low-luminosity spiral 
galaxies. 

Configurations obtained for positive values of nonminimal coupling exhibit positive pressures and hence are 
quite similar to the ordinary boson stars. In Fig. [9] we plot the energy density, radial and transversal pressure 
for A0 = on the left panel and A^ — 100 on the right panel. For each A^ two configurations are presented, one 
for £„^ c and the other one for a nonminimal coupling that is much larger then the critical one. For each 
the effect of increasing £ is only to decrease mass and increase radius (thus decreasing the compactness) without 
any drastic changes in the behaviour of energy density and pressures. However, the profiles of the energy density 
and pressures qualitatively do change considerably for different A^. On the right panel of Fig. [5] the hump in the 
energy density occurs as £ increases. This hump is actually followed by a violation of the strong energy condition 
which is more significant for larger £. Hence the hump in the inset of Fig. [U] is more prominent. In order to justify 
this statement in Fig. [10] we plot the energy density and pressures for a configuration that strongly violates the 
SEC (in the region of negative transversal pressure - see inset of Fig. [TUj) . 

In this subsection, apart from qualitative behaviour of the energy density and pressures, one could also infer 
subtle relations among the total masses and radii. In particular, increasing the central field value is followed by 
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FIG. 8: The energy density, radial (solid curve in the inset) and transversal pressure (dashed curve in the inset) for 
£ = -0.7, A = 0, cr = 0.1 and ml = 0.9 m 2 P . 
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FIG. 9: The energy density (solid), the radial pressure (short-dashed) and the transversal pressure (long-dashed) for 
a) left panel: A^, = and {£,<r c } = {0.6,0.2635} and in the inset {6.4,0.0364} and b) right panel: \$ = 100 and 
{£> M = I 7 - 8 . 0.1194} and in the inset {12.8, 0.0845}. Also ml = m%. 



a decreasing radius up to the maximally allowed mass. This interplay among the mass and the radius is best 
explored by analyzing the effective compactness. 



C. Effective compactness 

Following Ref. [48| . we define the effective compactness as 

C(o- ,A*) = M ™ i ° (hX *\ (37) 
-K99 

where Rgg is the radius at which the mass, defined in terms of the metric function e~ A = 1 — 2m/r, equals 99% of 
the total mass M — m(oo). The effective radius owes this sort of definition as the scalar field is (exponentially) 
infinitely extended and thus always with zero compactness. 

As shown in Ref. (48| the effective compactness in a minimal setting increases with the self-coupling and as 
— > 00 the maximal effective compactness approaches C max ~ 0.16 as shown in the right panel of Fig. Q] For 
each A<£, the maximal compactness corresponds to the parameters matching the critical field value o~ c . That is 
the maximally allowed mass and its radius. 

If M max is not constrained by the weak and dominant energy conditions, the effective compactness for £ > is 
largest for large £ and A^ = as shown in Fig. [11] and approaches C max ~ 0.20. This value is only slightly larger 
then the maximal effective compactness obtained in the minimal setting and can be related to a SEC violation. 
Why this value is not larger, probably can be explained with the fact that the SEC is violated only near the 
surface where the transversal pressures become negative. 
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FIG. 10: The energy density and the principal pressures (insets) for m| = m%, £ = 16, {\^,,a c } = {0,0.145}. 
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FIG. 11: The effective compactness as a function of the ^-coupling for A^, = (dotted curve), A^> = 20 (short dashed curve) 
and \<f, = 50 (long dashed curve). Also mi = m 2 P . 



However, for £ < the compactness is much greater and reaches its maximum value for large negative values 
of the nonminimal coupling and also in the case when A^ = 0, which approximately equals C max > 0.25. Even 
though the strong energy condition is not violated in this region, an increased effective compactness can be 
attributed to negative pressures that weaken gravity, thus enabling more matter to be accommodated in a fixed 
volume. This result is also very important as it suggests that, in order to build a highly compact object, we ought 
to have configurations with negative principal pressures and a violation of the SEC. 

When restrictions from the weak and dominant energy conditions are included, the effective compactness 



behaves as shown in Fig. [12j The maximal values for each A^ are obtained for £ crit 



WEC.DEC 



and then abruptly 



decrease with incresing/decreasing nonminimal coupling. This brings us to the conclusion that the most compact 
objects are produced in the domain of negative pressures and large self-couplings, with the maximum effective 
compactness only slightly larger then the minimal case C max > 0.16. 

Nevertheless, figures [Til and [T2l are very useful, as one can easily relate the mass to the radius. If we want to 
create a compact object with, for example, a radius of R — 15 km, then from C max one can easily read off its 
mass. The range of effective compactness C max = 0.05 — 0.25 correspond to the masses M = (0.5 — 2.5) Mq. 
However, in order to obtain scalar's masses and self-couplings, one needs to employ the universality described in 
Sec. II. By fixing the radius to, e.g., R — 15km one can calculate /3 for any configuration with differing radii 
(in the reduced Planck units). Then, the scalar's mass and self-coupling are easily obtained by applying the 
rescaling conditions m| — > m\/ ' (P and A^ — > A^//? 2 . By inspection of all above diagrams depicting the energy 
density and pressures (in previous subsection) it can be inferred that the radii of all given configurations roughly 
fall within the range r = (10 — 40) lp. Hence, if we want to create a star of R — 15 km the corresponding /3s 
are /3 = (18.6 — 4.6) x 10 36 leading to the scalar's masses = (0.27 — 1.08) x 10~ 8 eV, which could be in the 



13 




-5 5 10 



FIG. 12: The effective compactness as a function of £ for configurations that obey the WEC and DEC. A^ = (dotted 
curve), A^ = 20 (short dashed curve) and A^ = 50 (long dashed curve). Also m| = mp. 

range of the neutrino masses. To calculate the rescaled self-coupling, let us, for convenience take its starting value 
A^ = 50. After the rescaling we obtain A^ = (14 — 0.24) x 10~ 73 . But of course, if one considers a case when the 
coupling has reached saturation, one could increase the value of the un-rescaled A^> arbitrarily, which would then 
yield more reasonable (i.e. larger) values of the rescaled A^. 

V. CONCLUSIONS 

In this paper we have examined spherically symmetric configurations made of a scalar field nonminimally 
coupled to gravity. Our results are in perfect agreement with those of Colpi and collaborators (Ref. [l(|) for 
self-interacting minimally coupled scalars and with those of Bij and Gleiser (Ref. for a non-self-interacting, 
nonminimally coupled scalar field. 

We have showed that already a minimal extension of Einstein's theory to the nonminimal coupling results in 
radically different configurations from standard boson stars, i.e. dark energy- like stars which are characterized 
by negative principal pressures. Upon investigating the energy conditions in more detail, it turned out that 
the strong energy condition, which should be violated in the interior of dark energy stars, and whose violation 
signals repulsive gravity, is satisfied in the interior of these configurations. However, we presented an example 
of a proper dark energy star, i.e. with negative pressures and a violating SEC, which is obtained for a negative 
self-interaction. Even though the configuration presented here does not suffer from violation of the weak and 
dominant energy condition, configurations with negative potentials, in general, are not that appealing due to their 
stability issues. We also presented one higher mode solution that led to gravastar-like principal pressures. That 
is, both principal pressures reveal positive atmosphere (region near surface). But, the strong energy condition is 
obeyed while the weak and dominant conditions are violated, thus again without spacetime regions with repulsive 
gravity. 

When imposing restrictions on classical matter by energy conditions we found regions of parameter space for 
which both the weak and dominant energy conditions are violated. In particular, the weak energy condition is 
violated for all negative values of the nonminimal coupling, if it is less then a critical value £ < £St ■ The 
dominant energy condition is violated for all positive values of the nonminimal coupling if greater then a critical 
value £ > £„^ C - The consequences of a violation of the WEC and DEC are encoded in the maximally allowed 
masses that are now shifted to lower values. The strong energy condition is violated in the region of a positive 
nonminimal coupling and is followed by humps in the energy density. Even though violation of the energy 
conditions does not support the view of classical matter, it would be of interest to explore in more details the 
imprint of a test particle moving in such a background. 

Furthermore, we analyzed the effective compactness for configurations that do or do not satisfy the WEC and 
DEC, and found that the maximum effective compactness is attained in the regimes of negative pressures for non- 
self-interacting configurations and equals C max > 0.25 for configurations that violate the WEC and C max > 0.16 
for those configurations that obey the WEC and DEC. This result sets limits on the boson star mass. For example, 
when R = 15 km the maximum mass is M = (2 — 2.5) Mq which belongs in the domain of neutron stars. Even 
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though the strong energy condition is not violated, an increased maximum effective compactness could be related 
to the existence of negative pressures. 

In addition, we developed a universality condition based upon which one can calculate scalar's masses and self- 
couplings for all given configurations. Even though in this paper we focused on parameters that yield compact 
objects, with the universality condition it is possible to extend this analysis to larger structures that match 
galactic sizes, such as, for example, dark matter halos. An investigation of observational constraints on the model 
are underway. 

Although theories with a nonminimally coupled scalar field represent a simple and quite benign extension of 
general relativity, they provide a plethora of different interesting astrophysical structures, ranging from isotropic 
polytropes to highly anisotropic dark energy-like stars. Nevertheless, within this model it is, in fact, not possible to 
create a highly compact, nonsingular object whose characteristics are arbitrarily close to those of the Schwarzschild 
black hole. Yet, from this work one can infer that the real black hole mimicker might be produced in the context 
of modified theories of gravity. 
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